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Period-luminosity relations for Galactic Cepheid variables 
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ABSTRACT 

In this paper, we derive the period-luminosity (PL) relation for Galactic Cepheids with 
recent independent distance measurements from open cluster, Barnes-Evans surface 
brightness, interferometry and HST astrometry techniques. Our PL relation confirms 
the results of Tammann et al. (2003), which showed that the Galactic Cepheids follow 
a different PL relation to their LMC counterparts. Our results also show that the slope 
of the Galactic PL relation is inconsistent with the LMC slope with more than 95% 
confidence level. We apply this Galactic PL relation to find the distance to NGC 4258. 
Our result of = 29.49 ± 0.06(random error)mag. agrees at the ^ lAa level with 
the geometrical distance of /igeo = 29.28 ± 0.15rna(7. from water maser measurements. 
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1 INTRODUCTION 



Recently. iTammann et all ll2003l . hereafter T03) derived the 
Galactic PL relation by combining the Galactic Cepheids 
\ with independent distance m easurements from open clus- 
. ters/associations J^as^^l^ and from the Barnes- Evans 
(BE) surface brightness techniques iGieren et alll993) . The 
^ resulting Galactic PL relations in T03 are steeper than the 
■ LMC PL relations c ommonly applied in d istance scale ap- 
plications (as in, e.g.. lFree dman et alj|200 lll. S imilar conclu- 
sions are also reported in lFouau^^^i2003ll . 

The need to use the Galactic PL relation as a fundamen- 
tal cal ibrating relation has become mo r e desirable in recent 
vears JPeasd 120031: iFouaue et alJ 12003: iKanbur et alJl2003t 
iTammann et al.ll2003l: iThim et alJl2003ll . because of the foL 
lowing two main reasons: (a) The average value of metal- 
licity (defined as 12 -I- \og[0/ H]) in target galaxies of the 
Ho Key Project is ~ 8.84 ± 0.31de2; (Frccdman et al...200ll) . 
which is closer to the standard Solar value of 8.87± O.OTciea; 
JCrevesse et alJll996l) than the LMC va lue of 8.50±0.08dea; 
(see reference in iFerrarese et alll2000h : and (b) There is 
evidence that the LMC PL relation is broken at 10 days 
jTammann et al.l 120021: iKanbur fc Ngeowl liooi) . i.e., the 
short (< 10 days) and long period Cepheids in the LMC 
follow different PL relations. Due to these reasons, the cali- 
brated Galactic PL relation will become important in future 
distance scale studies. 

In this paper, we derive the Galactic PL relation from 
Cepheids with independent distance measurements. Our 
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analysis of the Galactic PL relation is similar to T03 but 
different in the following aspects: 



(i) In add i tion t o the Cepheids from iFeasd il999l) and 
ICieren et al.l il998fl that are used in T03, we include other 
recent distance measurements to Galactic Cepheids that are 
available in the literature (see Section 2). These include 11 
additional Cepheids that are not included in T03. 

(ii) As most of the Cepheids we considered here have more 
than one independent distance measurement, we took the 
standard weighted-mean of the available distances as the 
final adopted distance to derive the PL relation. 



2 GALACTIC CEPHEIDS WITH 
INDEPENDENT DISTANCE 
MEASUREMENT 

We collected the Galactic Cepheids with recent distance 
measurements from the literature, which include: 

(i) Distances from Open Cluster techniques: These 
distances are adopted from table 3 of T03, where the authors 
adopted a distan ce modulus for the Ple i ades of Upieiades = 
5.61 ± 0.03mag. dStello fc Nissenll200lD . iFeasd lll999D esti- 
mated that the uncertainty associated with cluster distanc e 
moduli is ~ 0.15-0.20mag. (see, e.g.. iRomeo et alJll989ll . 
hence we assign an uncertainty of Q.2Qmag. to the open 
cluster distance moduli {^o[O.C.]) in Table This error 
could incorporate the uncertainty due to the location of the 
Cepheid in the cluster (far-side vs. near-side), the uncer- 
tainty due to metallicity corrections from cluster to cluster. 
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etc. (see also iTurner fc Burkdiiooil . In addition, we also 
includ e oth er open clus t er dis tances from iTurner fc Burkd 
('2002') and ' Hovle et all l|200.tl in Table □ Since these ad- 
ditional open cluster distances are based on ^ipieiades ~ 
5.5(imag., we add a correction of +0.05mag to these dis- 
tances'^ in Table Q to be consistent with the T03 open clus- 
ter dist ances. Howev er, it is unclear whether the distances 
given in Feast ( 1999, where the T03 adopted distances origi- 
nate) and in Turner &J3urkc ( 2 0021) are totally independent 
to each other or not, since some of them were adopted from 
the same sources. For example, the open cluster distances 
for CV MON in Feast ( 1999) and Turner & Burke [20^ 
are all adopted fro m Turne r et al.l lll998lh . We have labelled 
the ITurner fc Burke (,2002i) distances that are dubious in 
this dependency regard in Table Q and excluded them in 
obtaining the weighted-mean distances. For GT CAR, CG 
CAS and TVCMa in Turner & Burke (2002), which are not 
included in lFeastI l)l999fl . there is no I band data available in 
the literature. We exclude them in order to have a consistent 
number of Cepheids in B, V and I bands. 

(ii) Distances from BE surface brightness tech- 
niques: We adopt these distances fr om Fou guc et al^ (20030, 
which is an updated version of FCieren et al] lll998l). with 
additional Cepheids that are not included in iGieren et alJ 
(^23)- Fouque (2003, private c ommunication ) has veri- 
fied that the the d istances from iGieren et alJ il998l) and 
iFouaue et alJ (l2003h are not independent of each other, 
as the distances from the later paper use the latest avai l- 
able data to update the distances in IGieren et al.l il99dl . 
The exception is CS VE L, where we include the distance 
from in ieren e t al. (199 ^ since there is no distance given to 
this Cepheid in Fouque et al. ll2003il. Not e that T03 only 

but not 



their paper. We did not include the distances from Hippar- 
cos because the errors in the distance moduli, after conver- 
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inclu ded the distances fr om 

from iFouaue et al.l i2003ft . We also include t he latest dis- 
tance measurements using BE techniques from lBarnes et all 
i2003h . who used a Bayesian statistical approach in their 
analysis. 

(iii) Distances from interfero metry tech n iques : 

These distances arc adopted from Lane e t alJ i2002f) . 
iNord grcn ct al. (2002) and Kcrvella et al. (200j). Note that 
for distances in lKervella et alJ ll|2003l) . we did not use the di- 
rect distance measurements for the four Cepheids {rj AQL, 
W SGR, l3 DOR & I CAR) because the errors are asymmet- 
ric and large (except for / CAR). Instead we adopt the dis- 
tances from their table 12. These are obtained from combin- 
ing the interfe rometry data and th e empirical Period- Radius 
relations from IGieren et all l)l99jj) . 

(iv) Distance from HST astrometry techniques: 
Currently there is only one Cepheid with a distance mea- 
surement fro m astrometry u s ing th e Hubble Space Telescope: 
S CEP from lBenedict et al.l (|2002h . 

The selected Galactic Cepheids and the corresponding 
distances from these sources are summarized in Table 
Since these distances are from independent measurements, 
we can take the weighted-mean among the available dis- 
tances (i.e. column 2, 3, 5 fc 6 in Table 0. Note that T03 
only uses the distance moduli in column 2 & 4 of Table Q in 



^ Thanks to G. A . Tammann for pointing out this correction to 
iHovIe et alJ i2003^ distances. 



sion from parallax, are large (see, e.g.. lMadore fc FreedmanI 
[l998). 

In this work, we exclude all the possible non- 
fundamental mode Cepheids and those not classified as 
"DCEP" in the Gener al Catalogue of Variable Stars (GCVS, 
iKholopov et alJll99^ . Some of them are mentioned in T03, 
and include: EV SCT, V1726 CYC, SZ TAU, QZ NOR, 
a UMi and V367 S CT. We further exclude EU TAU from 
iBarnes et alj ll2003h because it is a first-overtone Cepheid. 
Note that LS PUP is not classified as "DCEP" in the GCVS. 
We also exclude this Cepheid though it is incl u ded i n the 
T03 sample. For Ce pheids in ITurner fc Burkd tOQ'j) and 
iKervella et al.l (|2003f l , we did not include the Cepheids that 
are classified as "DCEPS" in the GCVS, which includes SU 
CAS, GH CAR and Y OPH. The additional Cepheids in- 
cluded in our sample but not in T03 are labelled in Table 

m 

In principle the number of Cepheids can be increased if 
we include the first overtone (FO) Galactic Cepheids in our 
sample, by using their fundamental mode periods in obtain- 
ing the PL relation. However, we prefer not to include the 
FO Cepheids due to the following reasons: (a) we want to 
avoid the contaminations from other types of Cepheids and 
select only the bona fide fundamental mode Cepheids; (b) 
we want to be consistent with T03, who excluded the non- 
fundamental mode Cepheids in their study; (c) the physics 
involved in fun damental mode and FO Ceph eids is n ot the 
same (see, e.g., lArrtmrfiUdLlQi^ lArrtmiel lo fc Aikaw^l^^j- 
iBono et a l.'l999a':' Feuchtinger et a l.'2000V' Bono et alJ2002^ 
and (d) the results of microlensing surveys to the Magel- 
lanic Clouds suggest th at the FO Cepheids follow their own 
PL relations (se e, e.g.. [ Udals ki et al."l999Y The theoreti- 
cal stu dies from lBono et air tl999b') and Baraffe fc Alibe^ 
i200l|l also suggest that the FO Cepheids follow difi'erent PL 
relations. Nevertheless, we list the FO Cepheids (classified 
as "DCEPS" in the GCVS) with recent distance measure- 
ments in Table 121 for completeness. Note that although BD 
CAS has been updated to "DCEPS" by Porctti (1994), there 
is no I band data available for this Cepheid in the literature. 
We therefore exclude BD CAS in Tabled 



3 THE PERIOD-LUMINOSITY RELATION 

The values of log(P), {B ~ V)o color, E{B - Vf, and the 
mean B, V and I band^ magnitudes for the Cepheids in 
Table0are all taken from T03. The absorption-to-reddening 
coefficient, R, for individual Cepheids is derived using the 
prescription given in T03^: Rv = 3.17(±0.13) + 0.44[(B - 
V)o - 0.78] + 0.05[E{B ~ V) - 0.42], Rb = Rv + 1.00 and 
Ri — Rv — 1.28. The B, V and I band extinction-corrected 



2 E{B - V)corr in T03. 

^ There are no mean I band values for rj AQL, S CEP and f 
GE M in T03, hence we adopted the I band mean magnitudes 
from lLanoix et al] 1199E|) . 

■* There are other formulae for the Ry available in the literature. 
We choose the formula of Rv from T03 in order to be consistent 
with the work of T03. The detailed study of the sensitivity of PL 
relation to the selected Rv will be addressed in a future paper. 
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Table 1. Distances to Galactic fundamental mode Cepheids. 



Cepheid 






/io(G'98)^ 




fj,o{other)'^ 


fio{w.m.)'' 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Tj AQL' 








6.986 ± 0.052 


7.526 ± 0.217''; 7.108 ± 0.148" 


7.025 ± 0.048 


RX AUR' 










11.101 ± 0.204'^ 


11.101 ± 0.204 


U CAR 


11.46 


11.46 ± 0.04S 


11.069 ± 0.038 


10.972 ± 0.032 




10.984 ± 0.032 


VY CAR 


11.63 


11.60 ± 0.09S 


11.419 ± 0.043 


11.501 ± 0.022 




11.503 ± 0.022 


WZ CAR 






12.980 ± 0.135 


12.918 ± 0.066 




12.918 ± 0.066 


I CAR 






8.941 ± 0.053 


8.989 ± 0.032 


8.670 ± 0.204«' 


8.981 ± 0.032 


CEa CAS 


12.69 


12.74 ± 0.15 






12.63 ± 0.14'i 


12.662 ± 0.091 


CEb CAS 


12.69 


12.74 ± 0.15 






12.63 ± 0.14'' 


12.662 ± 0.091 


CF CAS 


12.69 


12.74 ±0.15 






12.63 ± 0.14'' 


12.662 ± 0.091 


DL CAS 


11.22 


11.16 ± 0.04S 






10.99 ± O.U^ 


11.032 ± 0.115 


V CEN 


9.17 


9.16 ± 0.04S 


9.302 ± 0.024 


9.175 ± 0.063 




9.175 ± 0.060 


VW CEN 






13.014 ± 0.042 


12.803 ± 0.039 




12.803 ± 0.039 


XX CEN 






10.847 ± 0.065 


11.116 ± 0.023 




11.116 ± 0.023 


KN CEN 






12.911 ± 0.060 


13.124 ± 0.045 




13.124 ± 0.045 


5 CEP' 




6.76 ± 0.10 




7.084 ± 0.044 


7.173 ± 0.048''; 7.181 ± 0.089'' 


7.100 ± 0.029 


X CYG' 




10.30 ± 0.05 




10.421 ± 0.016 


10.209 ± 0.055<^ 


10.395 ± 0.015 


SU CYC' 




9.70 ± 0.06 








9.700 ± 0.060 


/3 DOR' 










7.566 ± 0.153° 


7.566 ± 0.153 


( GEM' 




7.85 ± 0.10 






7.794 ± 0.228''; 7.782 ± 0.211" 


7.732 ± 0.084 


Z LAC' 








11.637 ± 0.055 




11.637 ± 0.055 


T MON 


11.14 


11.05 ± 0.14S 


10.576 ± 0.067 


10.777 ± 0.053 


10.580 ± 0.068" 


10.721 ± 0.041 


CV MON 


11.22 


11.12 ± 0.04S 


10.901 ± 0.046 


10.988 ± 0.034 


11.39 ± 0.21'' 


11.003 ± 0.033 


UU MUS 






12.260 ± 0.092 


12.589 ± 0.084 




12.589 ± 0.084 


S NOR 


9.85 


9.82 ± 0.04S 


9.918 ± 0.025 


9.908 ± 0.032 




9.907 ±0.032 


U NOR 






10.769 ± 0.067 


10.716 ± 0.060 




10.716 ± 0.060 


TW NOR 


11.47 


11.47 ± 0.088 






11.38 ± 0.18'' 


11.393 ± 0.134 


V340 NOR 


11.17 


11.16 ± 0.11 


11.498 ± 0.130 


11.145 ± 0.185 


11.23 ± 0.12'' 


11.166 ± 0.070 


BE OPH 






9.496 ± 0.110 


9.271 ± 0.034 


9.265 ± 0.192'= 


9.271 ± 0.033 


UY PER 


11.78 


11.88 ± 0.50S 








11.780 ± 0.200 


RS PUP 


11.28 






11.622 ± 0.076 


11.160 ± 0.290'= 


11.555 ± 0.069 


VZ PUP 






13.551 ± 0.036 


13.083 ± 0.057 




13.083 ± 0.057 


AQ PUP 




11.78 ± 0.10 


12.750 ± 0.038 


12.522 ± 0.045 




12.397 ± 0.041 


BN PUP 






12.924 ± 0.051 


12.950 ± 0.050 




12.950 ± 0.050 


GY SGE 


12.65 


12.74 ± 0.088 


12.939 ± 0.071'' 






12.650 ±0.200 


U SCR 


9.07 


8.94 ± O.IOS 


8.869 ± 0.015 


8.871 ± 0.022 


9.137 ± 0.158'=; 9.13 ± 0.18'' 


8.881 ± 0.022 


W SGR' 










7.933 ± 0.169° 


7.933 ± 0.169 


X SCR' 










7.553 ± 0.161° 


7.553 ± 0.161 


WZ SGR 


11.26 


11.31 ± 0.04 


11.262 ± 0.021 


11.287 ± 0.047 


12.001 ± 0.169°; 11.23 ± 0.16'' 


11.316 ± 0.029 


BB SGR 


9.11 


9.08 ±0.08*5 


9.238 ±0.022 


9.519 ±0.028 


9.805 ± 0.181° 


9.518 ±0.027 


RY SCO 






10.469 ± 0.042 


10.516 ±0.034 


9.911 ± 0.147° 


10.485 ± 0.033 


KQ SCO 


12.36 


12.33 ±0.258 








12.360 ±0.200 


RU SCT 


11.60 


11.64 ±0.148 






11.41 ± 0.20'' 


11.480 ±0.141 


T VEL 






10.094 ±0.023 


9.802 ± 0.060 




9.802 ± 0.060 


RY VEL 






12.100 ±0.050 


12.019 ±0.032 




12.019 ±0.032 


RZ VEL 


11.19 


11.27 ±0.318 


11.169 ±0.025 


11.020 ±0.029 




11.024 ±0.029 


SW VEL 


12.08 


12.04 ± 0.058 


11.989 ±0.056 


11.998 ±0.025 




11.999 ±0.025 


CS VEL 


12.59 


12.59 ±0.148 


12.713 ±0.144 






12.671 ±0.117 


S VUL 


13.24 


13.24 ± 0.098 


13.731 ± 0.095'' 






13.240 ±0.200 


T VUL' 










8.920 ±0.146° 


8.920 ±0.146 


SV VUL 


11.83 


11.78 ±0.05 


12.325 ± 0.072'' 




11.331 ± 0.081°; 10.98 ± 0.21'' 


11.636 ±0.041 



fj,o(O.C.) = open cluster distance from T03; ^io{TB) = ope n cluster distance f rom lTurner fc Burkj j20o3), ad justed to fJ.pi,, i „dea = 
5.61mag. by adding +0.05mag.; fio{G98) = BE distance from lGieren et all il998l) : fio{F03) = BE distance from lFouQue et alj 12003) : 
lio{other) = distance from other sources; fio{w.m.) = the weighted-mean distance for the entries in column 2, 3, 5 & 6, when available. 
Excepts for CS VEL, which we include the fio{G98) in obtaining the weighted- mean. 

'' Distance measurements from interferometry; ij AQL &: ( GEM are from lLane et al.) i2002l) . and 5 CEP is from lNordgren et alj j20o2). 

° BE distance measurements from lBarnes et ahl l2003ll . who use a Bayesian approach in their analysis. 

'' Open cluster distances from lHovle et ahl l2003l) . adjusted to iJ-pidades = 5.61mag. by adding +0.05mag. 

° Interferometry distance measurements from iKervella et al.l i2003l) . 

^ iJSr astrometric measurement from lBenedict et al.l ^2002^ . 

8 These distances are not included in calculating the weighted-mean distances, see text for details. 

'' These distance moduli are not used in both lCieren et al] Jl99Si) and T03 as they appear to be outliers in the PL plots. 

i rrii — u-,;j„ 4- ; 1„J — 1 T^no 
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Table 2. Distances to Galactic first overtone Cepheids. 



Cepheid 


l^o(O.C.Y 




Mo(G98)=' 


Ato(S03)=' 


Mo(J^03)=' 


Mo(i^03)=' 


/.to(to.?n.)=' 


GH CAR 




10.99 ±0.20 










10.990 ± 0.200 


SU CAS 


7.12" 


7.11 ±0.03'' 










7.120 ±0.200 


V1726 CYC 


11.02 


11.02 ±0.03'' 










11.020 ±0.200 


QZ NOR 


11.17 


11.16 ± 0.11 


11.095 ±0.031 




11.23 ± 0.12 




11.169 ±0.075 


Y OPH 












9.06 ±0.21 


9.060 ±0.210 


EV SCT 


10.92 


11.09 ± 0.07 


11.066 ±0.033 




10.84 ±0.15 




11.026 ± 0.060 


SZ TAU 


8.72 


8.71 ±0.02^ 


8.090 ±0.042 


8.74 ±0.33 






8.725 ±0.171 


EU TAU 








10.27 ±0.16 






10.270 ± 0.160 


a UMi 


5.19 


5.44 ± 0.05 










5.425 ± 0.049 



^ fio{O.C.) = open cluster distance from T03; fioiT B) = open cluster fr om lTurner &: Burkd (2003), adj usted to Mp;p. i n,rfK« = 5.61mag. 
by adding +0.05mag.; ^o(G9 8) = BE distance f rom lFbugui et alj |20o3); fio(B03) = BE distance from lBarnes et alJ J^O^); ^10(^03) 
= open clust er distance from jHovle et all j20o3), adjusted to fJ.piciades = 5.61mag. by adding +0.05mag.; /io(^S'03) = interferometry 
distance from lKervella et alj l2003^ : Un(w.m.) = weighted-mean distances. 

'' These distances are not included in calculating the weighted-mean distances, since it is unclear whether these distances are independent 
to the distances given in Feast 1 1999) or T03. 

This distance is adopted from iPeastj (SS^), adjusting to f^pieiades = 5.61mag. by adding +0.04mag. 



absolute magnitudes for the 50 Cepheids in Tabled can be 
calculated by adopting the weighted-mean distance moduli, 
as given in the last column of Tabled and then fitted with 
least-square regressions to obtain the PL relation. The plots 
of the fitted Galactic PL relations are presented as solid lines 
in Figure with the following expressions: 

Mb = -2.594(±0.106) log(P) - 0.674(±0.123), a = 0.248 (1) 
Mv = -2.999(±0.097) log(P) - 0.995(±0.112), a = 0.226 (2) 
Mi = -3.303(±0.094) log(P) - 1.450(±0.108), a = 0.219 (3) 

The error bars in Figure are obtained from the 
quadrature sum of the error estimates of distance modulus 
(give n in the last colum n of Table Q, extinction (adopted 
from lFernie et all lll995ll database^) and mean magnitudes. 
We assign a conservative error of Q.05mag. to the mean 
magnitudes. This is reasonable because the mean magni- 
tudes are derived from accurate and reliable light curves 
jBerdnikov et alJl200(i) . 



1 1 1 1 1 
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3.1 The effect of FO Cepheids 

In this section we investigate the effect of including the 
Galactic PL relations if we include the FO Cepheids that 
are listed in Table |5| as suggested by anonymous ref- 
eree. The fundamental mode periods and the mean mag- 
nitud es in B, V and I band for these Cepheids are taken 
from iBerdnikov et all (|200(J) . The color excess for these 
Cepheids are calculated via the prescription given in T03: 
E(B - V) = 0.951g(-B - V)f, where E{B - V)f is taken 
from lFernie et alJ ^9^. Then the color, [B - V)o, can be 
calculated via {B — V) — E{B — V). The extinction correc- 
tions for these Cepheids are handled in the same way as in 
fundamental mode Cepheids. The fitted PL relations from 
the combination of 59 fundamental mode and first overtone 
Cepheids are: 



|http : / /ddo . astro ■ utoronto ■ ca/cephelds ■ html 



Figure 1. Galactic PL relation in B (top panel), V (middle panel) 
and I (bottom panel) band. The lines are the fitted PL relations, 
as given in equation (l)-(3). The solid circles are for fundamental 
mode Cepheids as listed in Table Q The open circles are the first 
overtone Cepheids from Table 1^ plotted with their fundamental 
mode periods. The error bars include the errors in distance mod- 
ulus, errors in extinction and an error estimation of 0.05mag. in 
mean magnitudes. 



Mb = -2.590(±0.100) log(P) - 0.665(±0.111), a = 0.265 
Mv = -3.007(±0.097) log(P) - 0.962(±0.108), a = 0.257 
Mi = -3.319(±0.098) log(P) - 1.406(±0.108), a = 0.260 

By comparing these PL relations to equations (l)-(3), 
we see that including the FO Cepheids does not significantly 
alter or improve the PL relations. However, the dispersions 
of the PL relations (a) have become larger than the PL 
relations without FO Cepheids. This is mainly due to the 
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one outlier, GH CAR (with log(Po) ~ 0.91), as shown in 
Figure After removing this Cepheid, the dispersions of 
the PL relations become comparable to those given in equa- 
tions (l)-(3). Even though the PL relations with the FO 
Cepheids are almost identical to the PL relations without 
FO Cepheids, we prefer the solutions given in equations (1)- 
(3) to be the calibrated Galactic PL relations, as we have 
argued in Section 2. 

3.2 The effect of open cluster distances from 
Turner & Burke 2002 

Since it i s unclear whether some of the open cluster distances 
given in iTurner fc Bur k^ l|2002|l are totally independent of 
the distances given in Ti^as? (Il999l) or not, we examine the 
changes of PL rela tions if we either exclude all open clus- 
ter distances from iTurner fc Burkj ll20nd) or assume that 
these distances are totally independent of lPeastl 1^9^, and 
include them in obtaining the weighted-mean distances. Re- 
call th at equations (l)-(3) used some of the Turner & Burke 
i2002r d i stance s that are either excluded in or independent 
of iFeastJ l^^^)- For the former case, the PL relations with 
49 fundamental mode Cepheids are; 

Mb = -2.627(±0.105) log(P) - 0.623(±0.123), a = 0.240 
Mv = -3.025(±0.096) log(P) - 0.954(±0.112), a = 0.220 
Mi = -3.320(±0.095) log(P) - 1.420(±0.110), a = 0.216 

For the lat ter case where we include all the open cluster 
distances from ITurner fc Burkj f2002), the fitted PL rela- 
tions for the 50 Cepheids in Table Q become: 

Mb = -2.619(±0.107) log(P) - 0.656(±0.123), a = 0.249 
Mv = -3.024(±0.097) log(P) - 0.977(±0.111), a = 0.225 
Mi = -3.328(±0.098) log(P) - 1.432(±0.108), a = 0.217 

T herefore, the exclusion or inclusion of lTxirner fc Burk3 
i2002h distances produces almost identical PL relations. 
These PL relations also agree and are consistent with equa- 
tions (l)-(3). However, there are som e open cluster distances 
that are not included in iFeasd \199^ or T03 (e.g. for Cepheid 
AQ PUP), and there are certain cases whe re we have confi- 
dence to b elieve that th e distances given in ITurner fc Bur'O 
(2002) and'Feasd il999ll are independent (e.g. for Cepheids 
CEa, CEb and CF CAS). Hence we continue to adopt equa- 
tions (l)-(3) as the calibrated Galactic PL relations. 



4 COMPARISON TO PUBLISHED RESULTS 

4.1 Comparison to other Galactic PL relations 

We selected recent empirical PL relations from the literature 
that give both V and I band PL relations to be compared 
with our results, because the PL relations in these two bands 
are e xten sively used in extragalactic distance scale studies 
(e.g., iFreedman e t aL , 2001). For exa m ple, we exclude the 
PL relations from lBarnes et alj J2003h or lHovle et ail ll2003ll 
because they did not give the I band PL relations in their pa- 
pers. The selected PL relations, along with our results, are 
given in Table |21 These include the Galactic PL relations 



derive d by ICieren et all (Il998l GAL-G981. iFouaue et alJ 
J2003l GAL-F03) and iTammann et alJ J2003l. GAL-T03). 
We exclude the Hipparcos-hased Galactic PL relations be- 
cause these PL relations adopted the slopes from LMC PL 
relations, and calibrate d the zero-points w ith Hipparcos data 
(as those used in, e.g.. iLanoix et al ] |l999l and lPaturel et alJ 
[2002). From the table, it can be seen that our results are 
consistent with the GAL-G98 and GAL-F03 PL relation. 
However there is some discrepancy between our results and 
GAL-T03. In this situation, we can us e the f-statistical test 
(see, e.g.. IZwillinger fc KokoskalbOOdl to assess the signif- 
icance of the difference in slopes under the null hypothe- 
sis that the slopes are the same. The results show that the 
p- value for the V and I band slopes are 0.15 and 0.27, re- 
spectively. Hence the null hypothesis cannot be ruled out at 
the 95% confidence level, and our results are also consistent 
with the GAL-T03 PL relation. 

The small discrepancy of the PL slopes between our re- 
sults and T03 is mainly due to the inclusion of additional 
distance measurements in this work (see Section 2 for de- 
tails). If we take the arithmetic, unweighted-mean of the 
39 Cepheids with fio{O.C.) and ^o(G98) in Table □ (as 
these distance moduli are used in T03), the slopes of the 
fitted PL relations become steeper: as ~ —2.750 ± 0.123, 
av = -3.130 ± 0.110 and ai = -3.402 ± 0.106, which agree 
with the results of T03. However, by taking the weighted- 
mean distances for these 39 Cepheids, the PL slopes become 
shallower but still consistent with T03: a_g = — 2.721±0.125, 
av = -3.102±0.112 and ai = -3.373±0.109. Thus the dif- 
ferent between our results and T03 is due to the inclusion 
of additional distance measurements in our study. 

For completeness, we also in clude the recent theoreti- 
cal PL relations in Table El from iBarafle fc Aliberti (1200 ll. 
GAL-Theoryl) and iFiorentino et alj ll200l GAL-Theorv2 
& GAL-Theory3) by adopting Z = 0.02. From the table, 
our empirical PL relations fairly agree wit h the theoretical 
PL relations from lBaraffe fc Aliberd i200 j) . However, none 
of the empirical PL relations given in Table|H]agree with the 
GAL-Theory2 and GAL-Theory3. 



4.2 Comparison to LMC PL relation 

Since the PL relation is shown to be different in the 
LMC and Galaxy by T03 (also in iFouaue et al.ll2003l and 
iKanbur et al]|2003l) . we verify this result by comparing our 
Galactic PL relation t o the LMC counterpart, as given in 
IFreedman et alj ||200 j) . For our Galactic PL relation, the 
difference in the V and I band slopes is: Aav = 0.239±0.102 
and Aai = 0.341 ± 0.096, which are 2.3o- and 3.6a results, 
respectively. We also apply the f-statistical test to test for 
the equality in the slopes of the Galactic and LMC PL re- 
lation. The results show that the slopes in the Galactic and 
LMC PL relation are inconsistent at more than a 95% con- 
fidence level, with p-value of 0.017 and 0.001 for the V and I 
band slopes, respectively. Therefore, Cepheids in the Galaxy 
and the LMC do follow different PL relations, and hence the 
Cepheid PL relation is not universal. 
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Table 3. Comparison of various PL relations'*. 



PL relation 


N 


ay 


bv 






bi 




Ref.*^ 


GAL-G98 


28 


-3.037 ±0.138 


-1.021 ±0.040 


0.209 


-3.329 ±0.132 


-1.435 ±0.037 


0.194 


1 


GAL-F03 


32 


-3.06 ±0.11 


-0.989 ±0.034 




-3.24 ±0.11 


-1.550 ±0.034 




2 


GAL-T03 


53 


-3.141 ±0.100 


-0.826 ±0.119 


0.24 


-3.408 ±0.095 


-1.325 ±0.114 


0.23 


3 


GAL-Here 


50 


-2.999 ±0.097 


-0.995 ±0.112 


0.226 


-3.303 ± 0.094 


-1.450 ±0.108 


0.219 


4 


GAL-Theoryl 




-2.905 


-1.183 




-3.102 


-1.805 




5 


GAL-Theory2 




-2.45 ±0.02 


-1.50 ±0.02 


0.17 


-2.78 ±0.01 


-2.02 ±0.01 


0.13 


6 


GAL-Theory3 




-2.22 ±0.01 


-1.62 ±0.01 


0.14 


-2.58 ±0.01 


-2.14 ±0.01 


0.10 


7 


LMC 


~ 650 


-2.760 ±0.030 


-1.458 ±0.020 


0.160 


-2.962 ±0.020 


-1.942 ±0.010 


0.110 


8 



^ My J = ay I log(P) + by J, and ay j is the rms dispersion. For LMC PL relation, assume fiLMC = 18.50mag. 

^ Reference: [1] Giereii et al.1 jl993); [2] iFouaue et al.l <2003D: [31 iTammann et alJ ^2003^ : [4] Equations (2) & (3) from this work; [5] 
Baraffo fc Alibert "<200 J) , with Z = 0.02; [6] table 6 from 'Fiorentin o et al.l 120021) . with Y = 0.31 and Z = 0.02; [7] table 6 from 
Fiorentino et alJ j20o3r ^with Y = 0.28 and Z = 0.02; fSl iFreedman et al.l <200ll) . 



5 THE DISTANCE TO NGC 4258 

The Galactic PL relation presented in the previous section 
can be used to find the distance to NGC 4258, because the 
m etallicity in this gala xy is 8.85 ± O.OGdex (see reference 
in lNewman et a^]l2001^ . which is closer to the Solar value. 
Furthermore, there is an accurate geometrical distance mea- 
surement to NGC 4258 using the water mase r in the in- 
ner disk of this galaxy jHerrnstein et aflll999^ . The mea- 
sured geometrical distance is 7.2 ± Q.5Mpc, corresponding 
to Hgeo = 29.28 ± 0.15mag. 

There are 15 Cepheids disc overed with HST obser- 
vations by iNewman et alJ i200 J) , where we adopted the 
periods and the mean V and I band magnitudes (from 
ALLFRAME photometry) for these Cepheids. The distance 
modulus to N GC 4258 can be obtained using the prescrip- 
tion given by iKanbur et al.. (,2003J. We di d not include 
the m etallicity correction because it is small iKanbur et alJ 
|20o3)- The results are given in Table |1] using the differ- 
ent PL relations in Table |3 This shows that the distance 
to NGC 4258 is consistent with these different empirical 
Galactic PL relations. However, the distance modulus of 
29A9mag. (from the PL relations given in equations [2] & 
[3]), and the median distance modulus of 29A6mag. (from 
the four empirical Galactic PL relations) is still ~ 1.4a and 
~ 1.2a, respectively, away from the water maser distance'^. 
Furthermore, all the theoretical PL relations in Table|3give 
a further and inconsistent distance modulus as compared to 
the distance moduli o btained from empirical PL relations 
(see, however, |Caputo ct al. 2002 for a way to resolve this 
discrepancy). This could be due, in pa rt, to the small numbe r 
of Cepheids discovered in NGC 4258 llNewman et al.ll200ll) . 
The on-going Cycle 12 ifST observations of NGC 4258 (Pro- 
gram ID: 9810; P.I.: L. Greenhill) that proposed to discover 
~ 100 Cepheids in this galaxy would help to solve the dis- 
crepancy between the Cepheid distance and water maser 
distance. 



^ Using the mean magnitudes derived in lKanbur et al.l JgOO^) for 
these 15 Cepheids decreases the fio in Table BI bv ~ Q.OGmag., 
with a median distance modulus of 29.40mag. 



Table 4. Distances to NGC 4258 with different PL relations''. 



PL Relation 






l-lo 


GAL-G98 
GAL-F03 
GAL-T03 
GAL-Here 


29.80 ±0.07 
29.79 ± 0.07 
29.73 ±0.07 
29.72 ± 0.07 


29.65 ± 0.05 

29.66 ± 0.05 
29.63 ± 0.05 
29.63 ± 0.05 


29.43 ± 0.06 
29.46 ± 0.06 
29.50 ±0.06 
29.49 ±0.06 


GAL-Theoryl 
GAL-Theory2 
GAL-Theory3 


29.80 ±0.07 
29.57 ±0.06 
29.41 ± 0.06 


29.74 ± 0.05 
29.58 ± 0.04 
29.45 ± 0.04 


29.66 ±0.06 
29.57 ±0.06 
29.50 ±0.06 


LMC 


29.90 ±0.07 


29.71 ± 0.05 


29.44 ± 0.06*' 



^ The errors arc ran dom error only. The systematic error is about 
0.l5mag (Newman c t alj|200ll) . 

^ The metallicity correcti on using the LMC PL relation to this 
galax y is Sz = +0.07mag. iFreedman et al ]|200i iNewman et alj 
|2003). Then the metallicity-corrected distance modulus becomes 
29.51 ± 0.06mag. 



6 CONCLUSION 

By using the recent independent distance measurements to 
50 Galactic Cepheids, we derive a Galactic PL relation in the 
B, V and I bands. Our analysis differs from that in T03 due 
to the following aspects: (a) we include other recent indepen- 
dent distance measurements; and (b) we took the weighted- 
mean to the available distances. The results confirm that the 
Galac tic PL relations are steeper than the LM C counter- 
parts iFouaue et alJl2003l: iTammann et al.ll2003l) . However, 
the steepness of the Galactic PL relation is still inconclusive 
(either as steep as in T03 or as shallow as in this study) 
because of the small number (~ 50) of Cepheids in the sam- 
ple. Application of the Galactic PL relation from equation 
(2) & (3) to determine the distance to NGC 4258 shows 
that there is still a ~ 1.4a discrepancy in distance with the 
water maser measurements. The Galactic PL relation can 
be improved if there are more independent distance mea- 
surements to Galactic Cepheids in the future, such as the 
Cycle 12 HST observations of the nearby Cepheids with as- 
trometric measurements (10 Cepheids from Program 9879, 
with P.I.: G. F. Benedict, where 3 of them are not included 
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in Table [Benedict 2003, private communication]) or the 
future Space Interferometry Mission (SIM, launch in 2009)'^. 
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